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Abstract 

Niobium diselenide has long served as a prototype of two-dimensional charge ordering, believed 
to arise from an instability of the electronic structure analogous to the one-dimensional Peierls 
mechanism. Despite this, various anomalous properties have recently been identified experimentally 
which cannot be explained by Peierls-like weak-coupling theories. We consider instead a model 
with strong electron-phonon coupling, taking into account both the full momentum and orbital 
dependence of the coupling matrix elements. We show that both are necessary for a consistent 
description of the full range of experimental observations. We argue that NbSe2 is typical in this 
sense, and that any charge-ordered material in more than one dimension will generically be shaped 
by the momentum and orbital dependence of its electron-phonon coupling as well as its electronic 
structure. The consequences will be observable in many charge-ordered materials, including cuprate 
superconductors. 
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The Charge Density Wave (CDW) order in 2i7-NbSe2 has been surrounded by controversy 
since its first observation nearly 40 years ago [1]. The most basic unresolved question 
concerns the driving mechanism of the ordering transition. In quasi-ID materials, it is well 
known that charge order arises from a Peierls instability: any interaction between states at 
Ep results in a periodic modulation of both the electronic charge density and the atomic 
lattice positions, which lowers the free energy through the corresponding opening of a band 
gap [2], Such a mechanism is termed ‘weak coupling’ since an arbitrarily small interaction 
will drive the system to order at low enough temperatures. In higher-dimensional systems, 
a nesting of the Fermi surface is typically required for such an instability to be effective. In 
NbSe 2 , however, Fermi surface nesting does not occur [1, 3-5]. The lack of a lock-in transition 
and the coexistence of CDW order with superconductivity below 7.2 K nonetheless seem to 
suggest a weak-coupling origin of the charge order. As a result a number of alternative weak- 
coupling mechanisms have been proposed, based on nested saddle-points in the electronic 
dispersion [6], local field effects [5], or a combination of weak nesting with momentum- 
dependent electron-phonon coupling [3, 7]. 

There are a number of experimental observations, however, which are hard to reconcile 
with weak-coupling approaches in NbSe 2 - The first concerns the size of the electronic gap 
in the charge-ordered state. Kinks in the Density Of States (DOS) observed by planar 
tunnelling experiments have been interpreted to arise from a gap of A = 35meV [8, 9], 
while only much smaller gaps between 2 — 5 meV were seen by high-precision Angle Resolved 
Photoemission Spectroscopy (ARPES) experiments [4, 10], and older studies even reported 
no gap at all in ARPES and resistivity measurements [11, 12], Assuming a weakly-coupled 
driving mechanism for the charge order, the size of the gap is expected to be proportional 
to the transition temperature, A ~ 1.76 ksTc, in direct analogy to the BCS theory of 
superconductivity [13]. A gap size of 35meV would then be far too large to explain the 
observed Xcdw of 33 K [14]. The 5 meV gap consistent with ARPES data is of the expected 
size, but it occurs only in isolated points on the Fermi surface [4, 10], raising the question of 
how the correspondingly small gain in electronic energy can overcome the cost of introducing 
the periodic lattice distortions associated with CDW order. Additionally, while the gaps 
seen in ARPES measurements are connected by the CDW wave vector, other equally well- 
connected points support either only a much smaller gap, or no gap at all [4, 10]. It was 
recently suggested, based on Scanning Tunneling Spectroscopy (STS) measurements [15], 
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that a particle-hole asymmetric gap of 12meV, centred above Ep, exists in NbSe 2 - The 
offset in energy explains why ARPES only sees a smaller gap size, while the value of 12meV 
is consistent with the order of magnitude expected for a strongly-coupled CDW, which may 
have A ~ 4 kpTc in analogy to strong-coupling superconductors [16-18]. 

The idea that the CDW order in NbSe 2 is of the strong coupling variety, in which a sizable 
electron-phonon coupling is essential in causing an instability of the electronic structure, is 
also consistent with two further recent experimental observations. Firstly, it was observed 
in soft X-ray scattering experiments that at the transition temperature a broad momentum- 
range of phonons simultaneously softens towards zero energy [7, 19]. This is in stark contrast 
to the sharp Kohn anomaly expected in a Peierls-like transition [20], and may be understood 
to be an effect of the importance of the entropy carried by phonon fluctuations in the 
strongly-coupled scenario [21, 22], These phonon fluctuations are localized in real space, 
and therefore necessarily involve a broad range of momenta. In an electronically nested 
system, on the other hand, the nesting vector localizes any fluctuations to a single point in 
momentum space. The second experimental result arguing for a strongly-coupled transition 
is the observation of a reduction of the density of states persisting well above the transition 
temperature, in what has been termed a pseudogap regime [4, 23, 24], In strongly coupled 
CDW materials, a pseudogap typically appears through the presence of locally fluctuating 
charge order without long-range coherence [23]. These same charge fluctuations have recently 
been suggested to play a central role in the description of the pseudogap phase of cuprate 
superconductors [25-28]. The precise nature of the pseudogap in these materials, however, is 
difficult to establish owing to the presence of various nearby competing orders. NbSe 2 , being 
free of such complications, provides a test bed on which the formation and characteristics 
of a pseudogap phase clue to charge fluctuations can be studied directly. 

Here we present a theoretical analysis of the charge order in NbSe 2 based on a model 
of strong electron-phonon coupling. Although the transition is phonon driven we End that 
a full knowledge of the electron states scattered from and to, including both wave vector 
and orbital dependence, is crucial in explaining the experimental data. Despite a lack of 
nesting, the limited matching of states at Ep acts in unison with the momentum dependence 
of the electron-phonon coupling to select out the CDW ordering vector, while the orbital 
characters of the bands naturally explain why a CDW gap appears primarily in only one 
of the Fermi surface pockets. The resulting momentum dependence of the gap itself agrees 
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with ARPES experiments as well as with the observed particle-hole asymmetry in the DOS. 
Finally we find that the phonon modes are softened over a broad range of wave vectors, 
and that including the presence of strong phonon fluctuations leads to a suppression of 
the transition temperature, implying a range of intermediate temperatures dominated by 
incoherent, fluctuating, charge order, and a corresponding pseudogap. In light of these 
findings, and the generic nature of the electronic structure of NbSe 2 , we argue that the 
momentum and orbital dependence of the electron-phonon coupling may be expected to play 
a similarly central part in the description of other charge ordered materials, especially in the 
presence of relatively strong coupling as in the case of the layered high-T c superconductors. 


RESULTS 

The extent of electronic nesting 

Niobium diselcnide is a layered material, in which hexagonal layers of Niobium atoms 
are sandwiched between similar layers of Selenium atoms, displaced so that they he above 
and below half of the Nb interstitial locations, as shown in the inset of Fig. 1. Consecutive 
sandwich layers are displaced to have the complementary half of the interstices occupied, 
giving two formula units per unit cell. The Fermi surface consists of two concentric barrel- 
shaped pockets centred around both the T and the K points, as well as a very small pancake- 
shaped pocket surrounding T, which we ignore in the analysis below. The inner pockets at 
K are seen in ARPES to develop the largest CDW gaps [4, 10]. The CDW itself is of the 
3 Q type, with three equivalent superposed modulations at 120° with respect to one other. 

To describe the electronic states, we employ a Slater-Koster tight-binding fit to the band 
structure using a basis of Se-4p and Nb-4d orbitals [29, 30]. For the two main bands crossing 
Ep that give rise to the barrel-shaped pockets (which we will refer to as inner, developing the 
CDW instability, and outer), the fit is constrained by data from ARPES measurements [10], 
while the remaining 20 bands are fit to recent LDA calculations [5]. In agreement with earlier 
reports [30], we End that both the inner and outer bands at Ep are primarily composed of 
the d> 3 Z 2_ r 2 orbitals centred on the two Nb atoms within a unit cell, in the form of bonding 
(inner) and anti-bonding (outer) combinations (as shown in Fig. 1). 

Focussing on the two non-interacting bands at Ep , we introduce the interaction between 
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FIG. 1. The band structure of NbSe 2 modelled by a tight-binding fit to both ARPES data and 
LDA calculations. The two bands crossing Ep (thick lines) are the ones involved in the formation of 
charge order. They consist of bonding and anti-bonding combinations of the two Niobium d% z 2_ r 2 
orbitals within a unit cell, as indicated. Upper inset: the Fermi surface resulting from the two 
bands crossing Ep. Lower inset: the layered atomic structure of NbSe 2 , with Nb in blue and Se in 
red. The unit cell, containing two formula units across consecutive sandwich layers, is indicated. 

electrons -0 and phonons (p: 


Hint = Ck+q fiqftk^k+q, (!) 

kq 

where the bare phonon dispersion can be described by a Brillouin function tending to a 
maximum of 11.2meV at the zone boundary [19]. The electron-phonon coupling g depends 
on the momenta, as well as the band indices m and n, of the the ingoing and outgoing electron 
states. Following Varrna et al. we model both of these dependencies using knowledge purely 
of the electronic structure, which has been shown to work well for a range of transition 
metals and their compounds [31]. In terms of the matrices of tight-binding overlap integrals 
5k and generalized eigenvectors the electron-phonon coupling matrix elements are given 


5 





































by [31]: 


—run _ m 

Skk' — v k 


'4-suk' 


Aj c S'k'^4k' 


' k') 


( 2 ) 


where a constant prefactor has been omitted, and the band velocity v n for electrons in band 
n is defined in terms of the band energy ££ as v£ = <9£k/<9k. In the interaction of equation (1) 
we restrict attention to the longitudinal component of g, since only the longitudinal acoustic 
phonons are observed to soften in inelastic neutron and X-ray scattering experiments [19]. 

Inserting the tight binding results into equation (2), we find that the inter-band electron- 
phonon coupling is strictly zero, while the intra-band coupling is about a factor of three 
stronger for the inner band than for the outer. This result immediately explains why the ex¬ 
perimentally observed CDW gap is so much more pronounced in the inner band [4, 10], since 
the static electronic susceptibility scales with the square of the electron-phonon coupling: 

_ f(£ m ) — f(f n ) 

V (q) « X Ig5+q| 2 p, _ (3) 

k ‘-’k+q Sk 

Here / ( E ) is the Fermi-Dirac distribution function, and the static Lindhard function is 
given by X 0 when g — 1. The absence of interband coupling causes the contributions of the 
two bands to the susceptibility to be entirely independent. This means that the lower value 
of the intra-band coupling in the outer band leads to an order of magnitude suppression of 
its contribution to X. 0j and therefore to a correspondingly reduced gap in its Fermi surface 
pocket. It is thus the orbital character of the electronic bands, and consequently the orbital 
dependence of the electron-phonon coupling, that causes a relative size difference between 
the CDW gaps in the pockets around the A'-point, even if they are equally well nested. 

The expression for the susceptibility in equation (3) can further be used to quantify both 
the amount of nesting present in NbSe 2 , and to evaluate the relative contribution of the 
electron-phonon coupling matrix elements in selecting the observed value of the (incommen¬ 
surate) ordering wave vector Q cdw = (1 — 5) § TM, with 5 ~ 0.02 [1], In the left panel 
of Fig. 2, three susceptibilities are compared. The dotted blue line shows the susceptibility 
for a perfectly nested model with a single ID band of cosine dispersion. The sharp peak 
characteristic of a nested Fermi surface is apparent. The dashed red line is the bare suscep¬ 
tibility in NbSe 2 , obtained by setting g = 1 in equation (3). This does not include any of the 
momentum and orbital dependence of the electron-phonon coupling, and is the form that 
has been widely used in previous studies [5, 10, 15]. It is almost completely flat except for a 
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FIG. 2. Left: electronic susceptibility as a function of momentum transfer. Blue dotted line: in a 
perfectly nested ID band structure, the susceptibility is sharply peaked at q = Q cdw- Red dashed 
curve: the result for the two bands in NbSe 2 crossing Ep, and an orbital-independent electron- 
phonon coupling. The susceptibility is nearly flat, without any peaks. Black solid curve: the 
susceptibility of just the inner band, in which a CDW gap is observed. Although a clear maximum 
can be seen at Q cdw-, it is not a sharp divergence indicative of perfect nesting. Right: evolution 
of the CDW vector Qcdw with varying peak position q 9 of the model electron-phonon coupling 
(see inset). Blue dotted line: in a perfectly nested ID band structure the electronic structure 
determines the CDW formation and Q cdw is independent of q 9 . Red dashed curve: the result for 
an orbital-independent electron-phonon coupling and the two bands crossing Ep, which is close to 
a straight line at 45°, indicative of a negligible role of the electronic structure. Black solid curve: 
the evolution relevant for NbSe 2 , given by just the inner band as imposed by the orbital-dependent 
electron-phonon coupling. This curve interpolates between the flat plateau of a perfectly nested 
structure, and the 45° line dominated by just the structure of the electron-phonon coupling. 

very broad and low maximum far from Q cdw-, an d can therefore not explain the observed 
ordering vector. The solid black line, on the other hand, displays the susceptibility of just 
the inner band. A broad but distinct maximum develops at the experimentally observed 
CDW vector. The selection of the ordering vector is thus shown to arise from the com¬ 
bined influence of both the electronic structure and the orbital and momentum dependent 
electron-phonon coupling, with neither contribution being negligible [4, 10]. 

Restricting attention to the inner band, we can quantify the extent to which nesting con¬ 
tributes to the selection of Q cdw- For this purpose we temporarily switch to a simplified 
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form of g, which is a function only of the momentum transfer: g = —a(q g — q ) 2 + g max - 
Even though in general the shape of the susceptibility is determined by both the incoming 
momentum and the momentum transfer, modelling the electron-phonon coupling in this 
simplified form allows us to directly quantify the effect of varying the position of its peak, 
and thus its influence on the position of the peak in Xq. The values of the parameters a, q g , 
and g m ax can be estimated by fitting the experimentally observed phonon softening [7] to the 
expression for the renormalized phonon energy in the Random Phase Approximation (RPA), 
^rpa = (^o — W>)> where we take fi 0 to be the high temperature (unrenormalized) disper¬ 

sion. In the right panel of Fig. 2 we plot the variation of the maximum in the susceptibility 
(i.e. the expected Q cdw) as q 9 is varied along TM). Again, the blue dotted line represents 
the situation of a perfectly nested band structure, in which the electron-phonon coupling 
never overcomes the peak in the susceptibility, and thus has no effect on Q cdw- The red 
dashed line has the same form of g for all bands, and is close to a straight line at 45° (grey 
dashed line). A true 45° line would indicate that Q cdw — q 9 , implying that the CDW 
vector is completely determined by the electron-phonon coupling, with the electronic struc¬ 
ture having no influence at all. The solid black line, which shows the behaviour of just the 
inner band, interpolates between these two extremes. There is no plateau, but there is an 
indication of some flattening, suggesting that while it is mainly the momentum dependence 
of the electron-phonon coupling which determines the value of CDW ordering vector, the 
role of the electronic structure cannot be neglected entirely. 


The CDW Gap 

Having identified the importance of both electron-phonon coupling and electronic struc¬ 
ture, we now revert to using the full expression of equation (2). To compare the size and 
shape of the CDW gap resulting from the orbital and momentum dependent g to the ex¬ 
perimentally observed behaviour, the RPA gap equations shown diagrammatically in Fig. 3 
are solved self-consistently. We simplify the calculation by using the 9x9 gap matrix 
which includes all higher harmonics appropriate for a commensurate CDW wave vector at 
Q = | T M, as an approximation for the gap matrix of the incommensurate CDW in NbSe 2 
with Qcdw — 0.98| TM. We further fix the values of the electronic self-energies, and con¬ 
strain the momentum dependence of the gap function to a five parameter tight-binding fit 



v (Go{k) if Qi = 0 (g(k) if Q t = 0 

k k+Qi |q if Qi 7^ 0 k k+Qi yT{k,k+Qi) if Qi ^ 0 



FIG. 3. Diagrammatic form of the self-consistent RPA gap equations. The single arrows indicate 
bare electronic Green’s functions Go, and the double arrows the renormalized Green’s functions Q 
and T. The anomalous functions E(k, k + Qi) = {^kWk+Q^ do not conserve crystal momentum, 
and correspond to the CDW order parameter. Notice that the final diagram includes a summation 
over the anomalous momentum Qi of the internal electronic propagator. 

consistent with the symmetries of the lattice. A self-consistent solution for the gap func¬ 
tion to all orders in RPA can then be found by searching for a fixed point in the flow of 
consecutive iterations of the gap equation. 

Fig. 4 shows energy and momentum cuts through the electronic spectral functions re¬ 
sulting from the solution to the gap equations, as well as their corresponding DOS. A good 
match with the shape of the experimentally determined DOS in STS experiments can be 
achieved by introducing a 4 meV shift in the chemical potential, which falls well within the 
range of experimental uncertainty [15]. The overall strength of the electron-phonon cou¬ 
pling, the only free parameter in the theory, is set to give a maximum CDW gap magnitude 
of « 12meV. The Fermi surface on the right of Fig. 4 then shows a CDW gap opening on 
the inner band only, and Fermi arcs forming at the same locations as observed in ARPES 
measurements [4, 10, 23]. The extent of the gap on the Fermi surface, its size (apparent 
from the dispersion’s back-folding in the centre of Fig. 4), and its restriction to one band, 
are all in agreement with experimental findings [4, 10, 23]. The plot of the DOS on the 
left indicates that these results, which were previously interpreted to imply a gap size of 
~ 5meV [4, 10], are in fact consistent with the STS results indicating a particle-hole asym¬ 
metric gap of ~ 12 meV [15, 23]. The asymmetry hides most of the gap structure above Ep, 
where photoemission intensity is suppressed by the Fermi-Dirac distribution [23]. The com¬ 
parison between the experimental STS data and the calculated RPA results shows that the 
particle-hole asymmetry arises naturally from the interplay between the electronic structure 


9 















FIG. 4. Left: the calculated density of states for the bare dispersion (black line) and renormalized 
dispersion (red line) compared to experimental results from STS measurements (blue points) [15]. 
The opening of the CDW gap leads to an asymmetric suppression of density of states, centred 
above Ep. Middle: the dispersion of the calculated spectral functions along MK. The back- 
bending induced by the gap at Ep is apparent in the band corresponding to the inner pocket 
about K. Right: Fermi energy cut of the spectral function resulting from the self-consistent gap 
equations, showing the formation of Fermi arcs as a result of gaps opening along MK. Below 
the dashed line the bare bands are shown for comparison. Note that the outer band develops no 
discernible gap; this can be compared directly to the results of recent ARPES measurements [4, 10]. 

and the momentum-dependent electron-phonon coupling. 


The Pseudogap 

The central role of the orbital and momentum dependent electron-phonon coupling in de¬ 
termining both the CDW propagation vector and the structure of the gap, suggests that the 
charge order in NbSe 2 falls in the strong-coupling regime, in which the electronic structure 
alone is not sufficient to characterize the emergent charge order. In the presence of strong 
electron-phonon coupling, it is well-known that the entropy associated with localized lattice 
vibrations plays an essential part in the system’s thermal evolution [21]. To understand the 
properties of NbSe 2 above its CDW transition temperature, it is therefore necessary to con¬ 
sider the influence of phonon fluctuations beyond RPA. The so-called Mode-Mode coupling 
Approximation (MMA) includes the lowest order terms of this type, and has been shown to 
successfully describe the influence of phonon fluctuations in some of the electronic properties 
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of related dichalcogenide materials like TaSe 2 and TiSe 2 [22, 32], The additional terms in¬ 
cluded in the MMA fall into two categories: a vertex correction, which we neglect by appeal 
to Migdal’s theorem because of the large momentum transfer associated with Q cdw, and 
self-energy corrections to the internal electron lines (indicated diagrammatically in Fig. 5), 
which we include. In agreement with recent X-ray scattering observations on NbSe 2 [7, 19], 
the MMA-renormalized phonon spectrum shown in the inset of Fig. 5 shows a broad range 
of phonon modes softening near the transition temperature, rather than the sharp Kohn 
anomaly indicative of well-nested materials. The momentum of the softest phonon mode at 
the transition temperature coincides with the experimentally observed CDW propagation 
vector. 

The softening of phonon modes, as the charge-ordered state is approached from high 
temperatures, is suppressed by the fluctuations beyond RPA. Since the temperature at 
which the first phonon mode reaches zero energy defines the transition temperature Tcdw, 
this means that the charge ordering transition temperature in MMA, Tmma, is suppressed 
with respect to the corresponding value in RPA, Trpa. Choosing the overall strength of 
the electron-phonon coupling such that T M ma matches the experimentally observed value 
of Tcdw = 33.5 K, Fig. 5 shows that the corresponding RPA transition occurs already at 
T rpa rj 300 K. Physically, this suppression of the transition temperature can be ascribed to 
the presence of fluctuations in the phonon field. In the temperature range Tmma < T < Trpa, 
the phonon fluctuations are strong enough to destroy the long range order predicted by (mean 
field) RPA theory. In that regime, the amplitude of the static order parameter vanishes, 
(Tq) = 0) but dynamic fluctuations of the order parameter persist, (|v?q| j ) > 0 [21, 22], 
The result is a locally-fluctuating short-range ordered state, characterized by the presence 
of a non-zero order parameter amplitude, without any long-range phase coherence [23]. Due 
to the non-zero amplitude, the gap in the electronic structure will survive even above the 
transition temperature Tcdw = Tmma, in the form of a pseudogap phase [4, 23]. The 
dynamical fluctuations in this regime may become locally pinned in the presence of defects, 
creating islands of order up to T = Tr PA ~ 300 K, in agreement with recent STM and X-ray 
scattering experiments [23, 24], The properties of this regime are reminiscent of the features 
which characterize the pseudogap phase of cuprate high-temperature superconductors [33], 
which have recently been shown typically to contain charge-ordered phases of their own [27, 
28]. Specific similarities are the presence of a pseudogap, Fermi surface arcs, dynamical 
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FIG. 5. Phonon dispersion in NbSe 2 - The solid black line indicates the bare phonon dispersion 
along TM, based on high temperature X-ray scattering results [19]. The blue dotted curve shows the 
dispersion resulting from the (mean field) RPA calculation at Tp pa = 303.5 K, while the red dashed 
curve represents the MMA prediction (including the phonon fluctuations shown in the diagram) at 
the same temperature. The fluctuations of the phonon field can be seen to suppress the onset of 
charge order predicted by the RPA theory. The bottom inset shows the MMA dispersion at Tmma = 
33.5 K, which represents the true transition temperature. It can be seen that the momentum of 
the softest phonon mode corresponds to the experimentally observed CDW propagation vector. 

fluctuations and a corresponding sensitivity to the timescales employed by experimental 
probes, as well as the existence of the crossover to the normal state at T* = Trpa- Note that 
the existence of a crossover temperature T* ~ 300 K, between the pseudogap and disordered 
regimes predicted here, falls outside of the temperature range that has been probed in NbSe 2 
to date [4, 24], 


DISCUSSION 

Although various theories for the formation of charge order in NbSe 2 have been able to 
explain many parts of its diverse experimental signatures, the material has so far evaded the 
simultaneous description of all of its properties by a single consistent theoretical account. 
The main features which remain difficult to reconcile with the weak-coupling theories consid¬ 
ered so far, include the emergence of CDW order in spite of a lack of nesting in the electronic 
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structure, and an anomalously flat electronic susceptibility, the softening of a broad range of 
phonon modes rather than a sharp Kohn anomaly, a discrepancy in the sizes of the CDW gap 
observed in different electron pockets in spite of them being similarly nested, disagreement 
about the existence, size and particle-hole symmetry of the CDW gap as seen by different 
experimental probes, and the presence of a pseudogap phase characterized by locally fluc¬ 
tuating charge order above the CDW transition temperature. We have shown here that all 
of these observations can be understood to result from the presence of a strong orbital and 
momentum dependent electron-phonon coupling. 

The combination of both the momentum-dependent coupling and the particular electronic 
structure of NbSe 2 is found to be crucial in the selection of the charge ordering vector, while 
the orbital dependence naturally explains the different behaviours of different Fermi surface 
pockets. The same dependencies on momentum and orbital character are also reflected in 
the emergence of a particle-hole asymmetric gap, which has a size of about 12 meV, centred 
above E F where it is largely inaccessible to observation by ARPES. These observations put 
NbSe 2 in the regime of strongly-coupled charge-ordered materials, in the sense that the 
specifics of the electron-phonon coupling qualitatively affect the properties of the ordered 
state. 

The pseudogap phase above T CDW may arise in such a scenario from local fluctuations 
of the lattice, which destroy the long-range phase coherence of the CDW order parame¬ 
ter but not its amplitude. Taking into account higher order phonon fluctuations, we find 
that indeed the mean-field transition temperature is suppressed, leading to a regime above 
Tcdw characterized by the absence of long-range charge order but the presence of phonon 
fluctuations. This agrees with the experimental observation of Fermi arcs and a pseudogap 
above the transition temperature, as well as that of pinned charge order at high tempera¬ 
tures around local defects. All of these signatures of the pseudogap phase are reminiscent of 
similar features seen in the pseudogap phase of cuprate high-temperature superconductors. 

In conclusion, while the Peierls mechanism can explain charge ordering in ID systems 
entirely in terms of the electronic structure, this is not the natural starting point for higher- 
dimensional materials. Instead, it is necessary to include full information about the de¬ 
pendence of the electron-phonon coupling on the momenta and orbital characters of the 
electronic states involved in the CDW formation. This is evidenced by the present case of 
NbSe 2 , a typical quasi-2D material without a strongly nested Fermi surface, in which knowl- 
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edge about all aspects of the electron-phonon coupling prove to be essential in explaining 
the full range of experimental observations. We suggest this consideration to be generically 
applicable across charge-ordered quasi-2D and quasi-3D materials lacking strong nesting 
features. In particular, the full momentum and orbital dependence of the electron-phonon 
coupling should be taken into account throughout the class of transition metal dichalco- 
genides and related material classes, including high-temperature superconductors. 
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